Medical Physics

Background. Accurate delineation of the gross tumour volume (GTV) - primary tumour (GTVp) and nodal disease (GTVn) - on FDG-PET/CT is a critical step of head and neck radiotherapy planning. Comparisons between lightweight custom networks and the auto-configured nnU-Net v2 are usually reported as end-to-end pipelines, conflating the contribution of the network with that of the inference-time post-processing applied on top of it. We separated the two. Methods. MiniUNet3D (custom 3D U-Net, 18.3 M parameters) and nnU-Net v2 (3d_fullres, 88.2 M parameters) were trained on the same 578 FDG-PET/CT cases (85/15 author-defined split of the HECKTOR 2025 Task 1 set, 8 centres) and evaluated on the same internal cohort. Three arms were compared pairwise: MiniUNet3D raw output at a fixed 0.5 threshold, MiniUNet3D with a locked adaptive post-processing pipeline, and nnU-Net v2. Comparisons used paired Wilcoxon tests with bootstrap confidence intervals, Bonferroni and Benjamini-Hochberg correction, and Cohen's d; catastrophic failure (Dice < 0.01) was compared with an exact McNemar test. Cases with an empty reference for a given target were excluded from that target's analysis (n = 98 GTVp, n = 93 GTVn). Results. With post-processing matched off, nnU-Net v2 was superior: median GTVp Dice 0.799 versus 0.592 (mean difference -0.244, 95 % CI -0.300 to -0.191; d = -0.88) and GTVn 0.774 versus 0.598 (d = -0.82). Post-processing raised MiniUNet3D to 0.800 (GTVp) and 0.738 (GTVn), recovering 79 % of that difference. Post-processed, MiniUNet3D matched nnU-Net v2 on GTVp Dice (p = 0.113) but remained inferior on nodal disease after Bonferroni correction (Dice p = 0.041; surface Dice p = 0.049). Catastrophic GTVp failures were 25/98 raw, 8/98 post-processed and 1/98 for nnU-Net v2 (McNemar p = 0.016). Inference took 34 s versus 78 s per case on the same GPU. Conclusions. Post-processing recovered most, but not all, of the difference between the two models, and it did not confer robustness: an eight-fold higher rate of empty contours on small primaries persisted, which is the more consequential difference for planning safety. Pipeline comparisons reported without a post-processing ablation risk attributing to a network what post-processing supplied.

Automated lesion segmentation is increasingly central to PSMA PET/CT interpretation, supporting staging, treatment planning, and response assessment at a scale that outpaces available nuclear-medicine expertise. However, automated PSMA-PET/CT whole-body lesion segmentation models are trained on images alone, with no knowledge of where in the body prostate metastases actually tend to occur. Radiologists use clinical domain knowledge of metastatic spread, but its absence in machine learning models produces false positives in anatomically implausible locations and missed lesions in high-risk sites such as the liver. In this work, we explore whether population-level spatial knowledge of metastatic spread can be used to augment deep learning segmentation predictions, and how such a prior should be fused with a network's output, without additional training. We build a data-driven metastasis atlas from 375 expert-annotated whole-body PSMA PET/CT scans and investigate its fusion with a trained segmentation network under a Bayesian framework, in which prediction probabilities from an nnU-Net-based lesion segmentation model serve as the likelihood and the data-driven atlas as the prior. Because metastases occupy only a small fraction of whole-body voxels, the atlas's peak probability is too low, and standard power-scaled or naive Bayesian pooling references lack the tools to deal with this shortcoming. This causes these standard fusion strategies to fail and, in the naive Bayesian case, to sharply degrade performance. We instead derive a calibrated, background-referenced log-odds fusion, one of many possible approaches to combine a population atlas with a deep learning model's predictions, distinct from classical multi-atlas label fusion in that it fuses a single population prior with a trained network's softmax rather than combining several registered atlases. Furthermore, this approach is neutral outside atlas support by construction, reduces exactly to the baseline network when unweighted, and requires no retraining. This atlas fusion significantly improved mean Dice over the baseline nnU-Net on a disjoint internal test set ($+0.011$, Holm-adjusted $p=0.021$) and on an independent external cohort ($+0.0129$, Holm-adjusted $p=3.8\times10^{-16}$), with lesion sensitivity improving from 0.849 to 0.861 internally and Dice improving over baseline in every stratified anatomic region, including the rare, high-risk sites motivating this work, while naive Bayesian pooling degrades performance sharply and power-scaled pooling underperforms it throughout. Our findings suggest that population-level spatial priors can meaningfully augment deep learning predictions in whole-body oncologic segmentation, provided the fusion rule is calibrated to where the prior actually carries signal.

Diffusion models achieve high reconstruction quality in low-dose computed tomography (LDCT), but their iterative sampling trajectories impose substantial computational costs. Unlike unconditional generation, paired LDCT reconstruction starts from an image that already contains the anatomy and spatial structure of the standard-dose CT (SDCT) target; reconstruction primarily requires correcting dose-related noise and artifacts. We therefore introduce Residual Endpoint Flow Matching (REFM), an LDCT reconstruction method that learns to transport an LDCT image directly toward its paired SDCT endpoint rather than defining a noise-to-image trajectory. REFM predicts the residual velocity along linear interpolations between both images and supports single-step and multi-step reconstruction using the same trained network. We evaluate five model capacities using 1 to 50 Euler steps against deterministic U-Net and diffusion-based baselines. Across all REFM variants, one-step inference consistently provides the highest reconstruction quality. On the TCIA validation set, REFM Base achieves 50.98 dB PSNR and 0.9865 SSIM at 94.54 fps, compared with 50.92 dB, 0.9847, and 9.26 fps for DDPM-10. REFM Small retains 50.71 dB while increasing throughput to 198.56 fps. Without fine-tuning, REFM Base also matches the 25-step DDPM baseline on the external Mayo Clinic dataset, although DDPM remains stronger on synthetically degraded CRLM images. Thus, our results show that exploiting paired anatomical correspondence enables diffusion-level LDCT reconstruction with a single step reconstruction.

Study design Retrospective cohort study. Objectives To correlate Hounsfield unit (HU) values, using elliptical regions of interest (ROI), that can be easily defined in routine clinical practice with magnetic resonance imaging (MRI) T2-hyperintense area fraction, as a surrogate for paraspinal muscle fat infiltration and to establish specific HU screening thresholds that may be applied with standard picture archiving and communication system (PACS). Methods We included 136 patients (71 men; 61.0 {+/-} 15.4 years) who underwent preoperative computed tomography (CT) and MRI within an 8-week period. Elliptical ROI HU values were measured at L2/3 and L4/5 for erector spinae, multifidus, and psoas major. MRI T2-hyperintense area fraction (Otsu thresholding) served as the fat infiltration reference. Linear mixed-effects (LME) models were used to assess the HU-T2 association and level-specific receiver operating characteristic (ROC) analyses (lower HU value side; n=136 per muscle-level) to identify thresholds for [&ge;]30% and [&ge;]50% infiltration criteria. Results Intraclass coefficients = 0.709 (HU) and 0.857 (T2 fraction); Goutallier weighted kappa = 0.579. In the overall LME, {beta} was -0.880 HU per 1% T2-fraction increase (95% confidence interval -0.935 to -0.825; marginal R2 =0.502); the association was steeper in multifidus ({beta} = -1.020) than in erector spinae ({beta} = -0.753). Psoas major (R = -0.226) was excluded from ROC analyses. Difference between L2/3 and L4/5 HU cutoffs was ~20 HU. The [&ge;]50% criterion revealed higher discrimination. Conclusions Elliptical ROI-based HU measurements may reliably screen paraspinal muscle fat infiltration in erector spinae and multifidus using standard PACS. Specific thresholds may allow practical preoperative evaluation without additional costs or radiation.

Dynamic positron emission tomography (PET) with long axial field-of-view (LAFOV) scanners enables multi-organ imaging and kinetic quantification beyond static (late-phase) imaging; however, the long times typically required for dynamic acquisitions remain clinically impractical. This study evaluates a deep learning (DL) framework to enable abbreviated dynamic PET acquisitions, comparing single-time-window (STW, early dynamic data only) and dual-time-window (DTW, early dynamic data plus a late 5-min static frame) protocols with early dynamic scan durations of 5-30 min and dose levels ranging from 360 MBq to 18 MBq. Seventeen 60-min dynamic [18F]FDG datasets were first motion-corrected using a staggered FALCON pipeline and then used to train and test a spatiotemporal DL model for autoregressive frame prediction. Performance was assessed across the full quantitative workflow, from DL-predicted frames and time-activity curves to organ-based kinetic modeling and voxel-wise parametric imaging in multiple tissues and two patient cohorts. DTW protocols consistently outperformed STW, better preserving late-phase kinetics. For a 15-min early dynamic scan, adding a late 5-min scan reduced mean absolute Ki difference from 23% (STW) to 17% (DTW) in the liver and from 26% to 15% in the thalamus. DTW + DL further reduced errors to [&le;]10% in the liver, thalamus, and breast lesion, and 16% in muscle. Our recommended protocol, 15-min early dynamic scan plus a 5-min late scan with DL, remained robust to up to a 5-fold dose reduction (~74 MBq). Overall, these findings support DL-enabled abbreviated, low-dose dynamic LAFOV PET as a clinically feasible approach for accurate kinetic quantification

Radiosensitizing nanoparticles represent a promising approach for enhancing the efficacy of proton radiotherapy; however, their performance is constrained by restricted penetration into tumor tissue, resulting in preferential perivascular accumulation. Here, we develop a spatially distributed mathematical model of a growing tumor undergoing proton therapy with intravenously administered radiosensitizing nanoparticles to investigate treatment optimization strategies. Using physiologically plausible parameter ranges informed by our own experimental measurements and published data, we demonstrate that co-administration of targeted nanoparticles with antibodies binding to the same tumor receptors can overcome transport-induced localization and promote a more uniform intratumoral redistribution of nanoparticles before irradiation. Population-level simulations across heterogeneous parameter sets suggest that moderate antibody doses consistently prolong tumor regrowth time, whereas higher antibody doses produce a pronounced and robust increase in tumor cure probability under a single high-dose irradiation regimen representative of preclinical settings. A key conceptual result of our analysis is the asymmetric risk associated with antibody co-administration. In contrast to antibody--drug conjugates, for which excessive dosing of unconjugated antibodies may severely compromise therapeutic efficacy, co-administration of antibodies with nanoparticle-based radiosensitizers constitutes a "safe-by-design" strategy with respect to tumor cell kill in the modeled single high-dose irradiation setting: although excessive antibody doses may yield suboptimal outcomes, they cannot reduce tumor cell kill below that achieved with targeted nanoparticles administered without antibodies. These findings identify antibody-mediated spatial redistribution of radiosensitizing nanoparticles as a favorable strategy that is expected to provide robust therapeutic benefit despite substantial variability in tumor characteristics.

Stroke remains a leading cause of mortality and long-term disability worldwide, yet rapid diagnosis is often limited by the shortage of trained radiologists, particularly in resource-constrained settings. Automated analysis of CT imaging offers a potential solution, but existing methods often struggle to achieve clinically generalisable performance while jointly addressing multiple diagnostic tasks. Here we present the Intelligent Integrated Stroke Diagnosis System IISDS, an end-to-end deep learning framework built upon StrokeGNN, a graph-based architecture that integrates 3D contextual feature extraction with U-Net-based 2D lesion segmentation to enable comprehensive stroke analysis from non-contrast CT scans. IISDS performs stroke subtype classification, lesion segmentation and lesion volume estimation within a unified pipeline. To develop and validate the system, we collected and curated BGD-ISD through a collaboration between AI researchers, neurologists, radiologists and clinicians, resulting in a large multi-centre dataset comprising 1,507 CT scans from 597 stroke cases acquired across six hospitals and medical centres in Bangladesh. Across BGD-ISD and multiple publicly available datasets, IISDS achieves state-of-the-art performance on all tasks, improving segmentation accuracy by [&ge;]0.011 Dice score, reducing lesion volume estimation error by [&ge;]0.3 average symmetric surface distance (ASSD), and increasing classification performance by [&ge;]0.018 area under the receiver operating characteristic curve (AUC) compared with existing approaches. These results demonstrate the potential of graph-based deep learning to enable clinically generalisable, automated and scalable stroke diagnosis from CT imaging, supporting rapid clinical decision-making, particularly in healthcare environments with limited access to expert radiological interpretation.

Despite the well-known health risks of neutron exposures, key gaps remain in understanding neutron-induced molecular responses and identifying reliable biodosimetric markers that distinguish neutrons from photon exposure. We provide the first genome-wide analysis of the human blood transcriptional response to an accelerator-derived fission-like spectrum of neutrons versus photons, evaluating transcriptomic relative biological effectiveness (RBE) and radiation quality-discriminating gene signatures. Whole blood from healthy donors was irradiated ex vivo with X-rays (140 kV, 0-4 Gy, n = 3) or neutrons (0.1-8 MeV, 0-1 Gy, n = 2), incubated for 6 h or 24 h, and processed for RNA sequencing from peripheral blood mononuclear cells (PBMCs). Neutrons were markedly more potent than X-rays at inducing differentially expressed genes (DEGs) at equal doses, showing a peak response 6 h post-irradiation followed by a decline. In contrast, X-rays caused a continuous increase in DEGs up to 24 h (neutrons vs. X-rays at 1 Gy: 1,449 vs. 121 DEGs at 6 h; 996 vs. 621 DEGs at 24 h). A universal p53-centered 34-gene signature, including FDXR, EDA2R, GADD45A, and ZMAT3, showed highly monotonic dose responses (Spearman correlation coefficient {approx} 1) across donors, radiation qualities, and timepoints. Additionally, difference-in-differences analysis identified radiation quality-discriminating genes only at 6 h, with transcriptional convergence observed by 24 h, suggesting a very narrow time window for biodosimetric differentiation. We identified a neutron-specific gene signature driven by cGAS-STING-NF-{kappa}B signaling (RELB, NFKB1, C3, MALAT1) and suppression of B-cell and myeloid identity genes (IGHD, TCL1A, CLEC7A, TLR2), defining a biologically coherent neutron quality index with distinct immunomodulatory effects. For the first time, we assessed neutron RBEs at the gene, pathway, and global transcriptomic levels in a human blood model, reporting a global transcriptomic neutron RBE of 1.30 (95% CI: 1.14-1.49) at 6 h and 1.21 (95% CI: 1.14-1.28) at 24 h, providing a valuable basis for biodosimetry in mixed-field exposure scenarios. Our findings advance the mechanistic understanding of neutron radiation responses and support the development of biodosimetric approaches for mixed-field exposure scenarios.

Purpose: The primary objective of the Stanford Knee Osteoarthritis PET/MRI Evaluation (SKOPE) study is to develop and evaluate a multimodal, dynamic [18F]NaF PET-MRI framework for characterizing whole-joint physiology and its relationship to osteoarthritis (OA) risk, pain, and disease progression. Specifically, we aim to integrate dynamic PET with quantitative and anatomical MRI, to characterize structural, compositional, and metabolic features across the knee and surrounding musculoskeletal system, evaluate acute tissue responses to exercise, and identify imaging biomarkers associated with OA risk, pain, and disease progression. Methods: The SKOPE study includes multimodal PET-MRI of the knee and surrounding musculoskeletal tissues, with imaging of the knee, tibia, ankle, thigh, hip, pelvis, and lumbosacral spine. Dynamic [18F]NaF PET is combined with conventional anatomical MRI and quantitative MRI techniques, including quantitative double-echo steady-state (qDESS) T2 mapping of cartilage, Dixon fat-fraction imaging, ultrashort echo time (UTE) T2* mapping of short-T2 tissues, UTE imaging of tibial bone, and zero echo time (ZTE) imaging for bone morphology and pseudo-CT generation. Additional MRI sequences characterize muscle composition, bone and joint anatomy, intervertebral discs, and regional vascular anatomy. Selected scans are acquired before and after a standardized exercise protocol to assess the acute physiological response of the joint. Automated segmentation is used to generate subject-specific masks of muscles, bones, vertebrae, and intervertebral discs. A subset of the MRI protocol is repeated at 1- and 2-year follow-up to assess longitudinal changes. Expected Impact: By combining dynamic bone metabolic imaging with quantitative measures of cartilage, menisci, muscle, bone, fat, vascular structures, and the spine and hip, the SKOPE protocol provides a whole-joint and multijoint framework for studying the structural, metabolic, and physiological processes associated with OA and pain. Exercise and longitudinal imaging further enable assessment of acute tissue responses and changes over time, supporting the development of quantitative imaging biomarkers for OA risk, pain, and disease progression.

Rib-cage morphology is a determinant of thoracic biomechanics, ventilation, and injury response, yet statistical shape models (SSMs) of the rib cage have relied on small cohorts (~100s of individuals) imaged by clinical computed tomography, which over-represents injury and disease. We constructed a surface-based SSM of the complete 24-rib cage from 26,275 standardised whole-body magnetic resonance imaging (MRI) scans of adults aged 19-74 years from the population-based German National Cohort (NAKO). Ribs were segmented with a deep-learning pipeline (a rib-extended SPINEPS model), reconstructed as per-rib surface meshes, and brought into dense vertex-wise correspondence by Gaussian-process morphable registration in Scalismo; the aligned ensemble was summarised by generalised Procrustes analysis and principal component analysis (PCA). Fourteen per-rib geometric descriptors provided a quantitative cross-walk between the abstract PCA modes and named shape features, and associations with sex, age, body size and composition (including body-fat percentage), and smoking exposure were estimated by multivariable regression with Benjamini-Hochberg false-discovery-rate control. Shape variation was strongly concentrated: 28 modes captured 95% of the total variance, and the first three alone accounted for 69.4% (PC1, 42.6%; PC2, 16.3%; PC3, 10.5%) and admitted consistent anatomical readings - a sexually dimorphic axis (PC1), a slender-versus-stout body-habitus contrast (PC2), and a free-rib-size axis at ribs 11-12 (PC3). The sexes were nearly fully separated along PC1 (Cohen's d = 2.52). Body mass and body-fat percentage were the dominant modifiable correlates of rib-cage shape, whereas the association with cumulative smoking exposure was comparatively small. The model is released as a population-representative geometric reference for benchmarking and morphing donor-derived finite-element human-body models and for further large-cohort shape analysis.

Objectives To validate two computational fluid dynamics (CFD) models derived from computed tomography angiography (CTA) for estimating trans-stenotic pressure gradients, using invasive intra-arterial pressure measurements as the reference standard in patients with renal artery stenosis (RAS). Background We assessed whether non-invasive assessment of the pressure gradient using CFD could be a reliable alternative to intra-arterial measurements for identifying hemodynamically significant RAS. Methods We performed intra-arterial measurements at rest and during dopamine-induced hyperemia to assess the trans-stenotic pressure gradient in 28 patients with RAS. A pre-intervention CTA scan was used to simulate the pressure gradient with a CFD model using a strategy based on Murray's law (CFD-Mu) and cortical volume (CFD-C). The agreement between the simulated and measured pressure gradients was assessed using intraclass correlation coefficients (ICC), Bland-Altman analysis and diagnostic agreement on the presence of a hemodynamically significant stenosis. Results In 20 patients, successful measurements and simulations were obtained. The ICC between measured pressure gradient and the CFD pressure gradient was 0.78 and 0.94 during baseline and 0.86 and 0.72 during hyperemia, for CFD-Mu and CFD-C, respectively. The sensitivity of CFD-Mu and CFD-C was 70% for both models at rest and 100% compared to the hyperemic measurements, whereas the specificity was 90% and 70% at rest and 79% and 72% during hyperemia, respectively. Conclusions The results support the use of individualized CFD simulations for hemodynamic assessment of RAS using CTA as input. The CFD models demonstrated high accuracy for the identification of a hemodynamically significant stenosis.

Ultrathin endoscopy is highly attractive for real-time tissue imaging in narrow and hard-to-reach regions of the body. A single multimode fibre (MMF) is an attractive probe because of its small diameter, flexibility, and diffraction-limited spatial resolution enabled by the large number of transverse modes guided within a single core. Because the distal fibre tip is inaccessible during endoscopy, reflection-mode imaging, in which the same fibre delivers illumination and collects backscattered light, is more practical than transmission-mode imaging. However, image recovery from the resulting speckle pattern is challenging because light undergoes double-pass propagation through the MMF, with mode coupling and dispersion; the backscattered signal is weak, and the camera records intensity only, without phase information. Here, we propose a single-shot reflection-mode MMF imaging framework that combines a reflected real-valued intensity transmission matrix (reflected-RVITM) with an image restoration network. The reflected-RVITM is calibrated using intensity-only measurements, without interferometry or phase retrieval, and provides a physics-guided initial reconstruction from a single backscattered speckle frame. A restoration network then refines this initial reconstruction instead of inverting the raw speckle. Four restoration backbones are evaluated: HPM-Attention-UNet, GAM, MambaIRv2, and CICPNet. On matched datasets, hybrid models outperformed corresponding networks trained to map raw speckle directly to images. For example, HPM-Attention-UNet on MNIST improved mean PCC from 0.572 to 0.944 (+65.1%). Under domain shift, with training only on Fashion-MNIST and tested on unseen CIFAR scenes, hybrid models achieved mean PCC of 0.61-0.65, compared with 0.36-0.50 for direct learning. This framework is further demonstrated using physical objects at the distal fibre tip. These results demonstrate that a reflected-RVITM physics prior combined with a restoration network enables single-shot image recovery after intensity-only calibration, offering a phase-retrieval-free and generalisable route towards minimally invasive reflection-mode MMF endoscopy.

Purpose/Objective There is paucity of data reporting outcomes of breast cancers with initial internal mammary nodal involvement and no visceral metastases, treated with curative hypofractionated radiotherapy . We report the outcomes from a tertiary centre alongside spatial patterns of recurrences in the above group Material/Methods For this retrospective cross-sectional study, consecutive patients contoured as per the ESTRO 2013 guidelines, treated between 2016-2022 were eligible if their diagnostic imaging demonstrated involvement of the internal mammary nodes. Radiotherapy (40 Gy/15#/3 weeks) was delivered to the residual breast / thoracic wall, SCF region corresponding to the ESTRO lymph node level 4 and internal mammary chain nodes. Residual IMN/ level 4 nodes received a boost of 10Gy/5#. Spatial mapping of sites of recurrence at the local site and three nodal sites (axilla, SCF and IMN) was performed using deformable image registration. Sites of recurrence at the local site and three nodal levels were contoured separately. Volumetric intersection of the recurrent gross tumour volume (GTV_recurrence) with treated clinical target volume (CTV) was calculated. Actuarial overall (OS), disease free survival (DFS) & cumulative incidence of local (LR), regional (RR) and loco-regional recurrence(LRR) were calculated using Kaplan Meier method. Univariate comparison of outcomes with or without residual disease was performed using the log rank test. Results The median age of the 61 eligible women was 49 years. 77% received neoadjuvant chemotherapy and the rest adjuvant chemotherapy. 82% patients had a mastectomy. Axillary lymph node dissection was done in 96.7%. Boosts to residual IMN and SCF nodes were delivered to 21(34.4%) and 2 (3.3%) respectively. Median follow up was 3.6 years. Out of the 61 patients, 42 patients were disease free with an estimated 3 year disease free survival of 75% (95% CI 64, 88%). Spatial mapping of locoregional recurrence was possible in all but 1 patient with local (only) recurrence who was lost to follow-up after mammogram only. Among the patients with loco regional recurrence 1 had recurrence in local site + SCF +axilla, 3 had recurrence in the SCF+axilla, 2 in the SCF+IMN and 1 in the axilla+SCF+IMN. Only one patient had isolated axillary recurrence or isolated SCF recurrence. There were no IMN only recurrences. Among the 8 patients with nodal recurrence, a total of 27 individual GTV_recurrence were identified in the axilla(n=11), SCF(n=11) and IMN (n=5). IMN recurrences showed complete or partial overlap with CTV. SCF recurrences were a mix with predominantly in-field recurrences while axillary recurrences occurred outside the treated volume.Four (6.6%) patients had Grade 2 lymphoedema as documented late side effect. Conclusion Aggressive treatment of IMN disease with adjuvant radiation is effective with good locoregional control. Systemic recurrences are common and may benefit from intensification strategies.

Peripheral artery disease (PAD) spans a continuum from large-vessel obstruction to distal microvascular dysfunction, yet routine non-invasive tests, including the ankle-brachial index (ABI), do not provide structurally resolved assessment of the foot microvascular bed and may be unreliable in the setting of medial arterial calcification or perioperative follow-up. Here we developed a clinic-oriented multispectral compound-scanning photoacoustic tomography system (MCPATS) for compression-free distal toe imaging, and an interpretable photoacoustic tomography distal microcirculation score, termed PACT-DMS, for phenotyping PAD-related distal vascular abnormalities. PACT-DMS was derived from anatomically standardized distal toe sections and integrated seven prespecified vascular features spanning trunk-vessel morphology, microvascular distribution and pulsation-related dynamics through a traceable linear support vector machine. In a prospective single-centre cohort of 45 participants, the bilateral fusion PACT-DMS model distinguished patients with PAD from healthy controls with an area under the receiver operating characteristic curve of 0.964 (95% CI, 0.907-1.000) and an accuracy of 91.1% (95% CI, 82.2%-97.8%) under subject-level leave-one-out cross-validation, supported by complementary robustness analyses. Exploratory analyses further showed that PACT-DMS identified abnormal distal vascular phenotypes in 6 of 9 clinically diagnosed PAD limbs with non-abnormal ABI and visualized distal vascular-bed changes before and after revascularization. These findings support MCPATS-enabled interpretable photoacoustic vascular phenotyping as a candidate adjunctive approach for distal microcirculatory assessment in PAD; larger multicentre studies with external validation and prespecified analysis protocols will be required to define its clinical role.

Rationale: Neuroblastoma is a devastating pediatric malignancy, for which surgical resection is a key factor in long-term survival. However, there are significant challenges in its resection, particularly in high-risk disease, as neuroblastoma encases surrounding critical structures, is often difficult to distinguish from desmoplastic or scar tissue, and can carry occult deposits of disease not readily identified on preoperative imaging or intraoperative visualization. Building on the principles of fluorescent and radio-guided surgery, in combination with the known overexpression of GD2 in neuroblastoma, we sought to develop and optimize 111In-Dinutuximab-IRDye800, a dual-modality GD2-targeted intraoperative molecular imaging agent, for use in pediatric neuroblastoma to help enhance patient safety while facilitating a more complete resection. Methods: Dinutuximab was conjugated to IRDye800 and DTPA, then radiolabeled with Indium-111 to yield 111In-Dinutuximab-IRDye800. Optimization occurred through ELISA assay to assess binding affinity, fluorescence intensity analysis to determine the optimal fluorescent degree of labeling, and phototoxicity testing through flow cytometry. Rodent models of neuroblastoma were then generated through injection of SK-N-BE(2) human neuroblastoma cells into the left adrenal glands of nude mice or RNU rats. A series of fluorescent and gamma biodistributions was performed, varying the dose, timing, and specific activity of the tracer. Tumor and organ uptake of the tracer was compared with one- or two-way ANOVA as appropriate, with Sidaks multiple comparison test to compare tumor uptake to individual organs. Once optimization was complete, a clinically significant events study modeled after human clinical trials was performed to evaluate the in vivo capabilities of 111In-Dinutuximab-IRDye800. Results: Increased ratios of IRDye800 per antibody led to decreased binding affinity for GD2 and was associated with formulation instability without significant return on fluorescence intensity. Specific activity of the tracer was not found to impact overall biodistribution of the tracer. A 45-50 microgram dose of 111In-Dinutuximab-IRDye800 with ratios around 1 DTPA and 1-1.5 IRDye800 per antibody imaged 4 days after tracer administration was found to be the optimal combination that maximized detectable tumor-specific signal. In the clinically significant events study mirroring human IMI clinical trials, fluorescent guidance identified additional malignant lesions not originally detected under white light in 64% of rodents. Conclusions: 111In-Dinutuximab-IRDye800 is a dual-modality GD2-targeted intraoperative imaging agent that is well-poised for clinical translation. As it preserves tumor specificity, yields clinically meaningful radiofluorescent signal, and is well-tolerated without adverse events after optimization was completed, it carries the potential to positively impact the safety and completeness of neuroblastoma resection.

Muscle architecture shapes muscle function and changes with age, growth, training and disease, yet quantifying three-dimensional (3D) muscle architecture in vivo remains challenging. We introduce a hybrid fascicle tractography approach for freehand 3D ultrasound data that accurately reconstructs 3D muscle fascicles with respect to an objective, anatomically relevant coordinate system defined by the muscle's central aponeurosis. The hybrid approach combines Hessian-based fascicle detection with wavelet-based refinement to generate volumetric fascicle orientations. In a synthetic dataset with known ground truth, fascicle orientations and lengths were estimated with errors of [&le;]2{degrees} and ~1.5%, respectively. In vivo, the approach detected physiologically plausible fascicle lengthening in the human tibialis anterior following a passive plantar flexion rotation, whereas diffusion tensor imaging of the same muscle did not. The proposed method enables anatomically relevant, objective and non-invasive quantification of 3D muscle architecture in vivo, providing a practical framework for applications in clinical and applied muscle physiology.

Background: Acute phase response is an early immunometabolic response to brain injuries, primarily coordinated by the liver via the activation of acute phase proteins. These immune responses can be both beneficial, promoting tissue repair, and detrimental, exacerbating neurological deficits, if not properly controlled. Despite the central role of the liver in immunometabolism, how hepatic metabolism dynamically adapts to traumatic brain injury remains under explored, primarily due to limited liver-specific modalities that can assess metabolic pathways in vivo. 13C MRI utilizing hyperpolarized 13C-pyruvate can assess key regulatory enzyme activities in hepatic metabolism. Methods: Rats with controlled cortical impact were studied in vivo using hyperpolarized [1-13C]pyruvate and [2-13C]pyruvate under fed and fasted conditions 3-4 days after injury. Hyperpolarized 13C products, including [13C]bicarbonate from [1-13C]pyruvate and [5-13C]glutamate, [1-13C]acetyl-L-carnitine, and [2-13C]phosphoenolpyruvate from [2-13C]pyruvate, were evaluated to assess mitochondrial and gluconeogenic metabolism. In parallel, liver tissues were collected following [U-13C3]pyruvate injection for NMR isotopomer analysis of phosphoenolpyruvate, glucose, and glutamate. Results: While no metabolic differences were detected under fed condition, [13C]bicarbonate and [2-13C]phosphoenolpyruvate increased after brain injury under fasted condition, indicating an upregulation of the hepatic gluconeogenic pathway after injury. 13C NMR of liver tissue extracts from injured rats showed an elevated [2,3-13C2]glutamate-to-[4,5-13C2]glutamate ratio and increased 13C-labeling in phosphoenolpyruvate than controls, confirming enhanced hepatic gluconeogenic pathway. Conclusion: This study demonstrates that hepatic acute phase response to brain injuries can be monitored in vivo by hyperpolarized pyruvate, which may be further utilized for longitudinal immunometabolic evaluation of the liver during pathogenesis and therapeutic interventions.

Patients with history of neck radiation for Human Papilloma Virus (HPV)-associated head and neck cancer experience rapid progression of carotid artery stenosis. The present study sought to determine whether screening and treating asymptomatic carotid artery stenosis in patients with a history of neck radiation is cost-effective. This study is a cost-utility analysis using a probabilistic Markov model over a thirty-year time horizon assessing carotid screening and treatment to avoid neurologic consequences of neck radiation for HPV-associated head and neck cancer. A strategy of no carotid surveillance was associated with a 14.8% cumulative risk of stroke and a strategy of ultrasound surveillance and treatment with TCAR was associated with a 3.0% cumulative risk of stroke. The latter had a median incremental cost of $1.04 million USD and provided a median 39.1 additional QALYs, resulting in a median incremental cost-effectiveness ratio of $26,556 per QALY. In conclusion, this study suggests that ultrasound surveillance and treatment with TCAR for asymptomatic carotid artery stenosis is likely to be cost-effective for patients who have been successfully treated with radiation therapy for HPV-associated head and neck cancer.

Diffusion-weighted MRI, beyond the commonly used diffusion tensor framework, offers a unique window into tissue microstructure in vivo, yet its clinical adoption has remained limited. Major barriers include the complexity of diffusion MRI sequence design, lengthy acquisition protocols, and the challenges associated with robust estimation of high-dimensional microstructural model parameters. Here, we address these limitations by combining optimised diffusion encoding with state-of-the-art simulation-based inference, establishing a clinically feasible framework for multi-compartment diffusion modelling. We validate the approach through i) in-depth in silico experiments and ii) in vivo studies made up of both human and rodent data. The resulting microstructural metrics are robust, reproducible across healthy individuals and show significant spatial associations with brain-wide expression patterns of cell-specific genes. Requiring less than 10 minutes of acquisition time, this framework substantially lowers the barriers to advanced microstructural imaging, a prerequisite step toward its eventual evaluation for the diagnosis, stratification, and monitoring of brain disorders.

Clinical brain MRIs from pediatric health systems represent a viable resource for modeling early neurodevelopmental trajectories and studying neurodevelopmental risk in real-world populations. However, a limitation to date has been the performance of existing segmentation tools for measuring various brain phenotypes in clinical scans. In particular, many tools underperform in infant scans due to morphological and physical changes such as rapid myelination. Here, we introduce ClinSeg: a robust segmentation approach tailored to early-life clinical MRIs with variable orientation, resolution, and contrast. We leverage existing registration and synthetic data generation tools to construct a training corpus for a 3d U-Net spanning anatomical and contrast diversity, including scans with morphological abnormalities from a pediatric hospital. Validated against manual segmentations, ClinSeg outperforms existing models in infancy while matching them in childhood and adolescence. Finally, ClinSeg enables the construction of reference brain growth trajectories in 11,699 individuals from 0-21 years of age, leading to the detection of more nuanced age-related findings in clinical groups.